arXiv:1504.07598vl [astro-ph.SR] 28 Apr 2015 


Astronomy & Astrophysics manuscript no. nissen.twins © ESO 2015 

April 29, 2015 


High-precision abundances of eiements in soiar twin stars ^ 

Trends with stellar age and elemental condensation temperature 

P. E. Nissen^ 

Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, 
Denmark, e-mail: pen@phys. au. dk. 

Received 7 April 2015 / Accepted. 


ABSTRACT 

Context. High-precision determinations of abundances of elements in the atmospheres of the Sun and solar twin stars indicate that 
the Sun has an unusual low ratio between refractory and volatile elements. This has led to the suggestion that the relation between 
abundance ratios, [X/Fe], and elemental condensation temperature, Tq, can be used as a signature of the existence of terrestrial planets 
around a star. 

Aims. HARPS spectra with S/A i 600 for 21 solar twin stars in the solar neighborhood and the Sun (observed via reflected light from 
asteroids) are used to determine very precise (cr 0.01 dex) differential abundances of elements in order to see how well [X/Fe] is 
correlated with Tq and other parameters such as stellar age. 

Methods. Abundances of C, O, Na, Mg, Al, Si, S, Ca, Ti, Cr, Fe, Ni, Zn, and Y are derived from equivalent widths of weak and 
medium-strong spectral lines using MARCS model atmospheres with parameters determined from the excitation and ionization bal¬ 
ance of Fe lines. Non-LTE effects are considered and taken into account for some of the elements. In addition, precise (cri 0.8 Gyr) 
stellar ages are obtained by interpolating between Yonsei-Yale isochrones in the log g-Tefs diagram. 

Results. It is confirmed that the ratio between refractory and volatile elements is lower in the Sun than in most of the solar twins (only 
one star has the same [X/Fej-Tc distribution as the Sun), but for many stars, the relation between [X/Fe] and Tq is not well defined. 
For several elements there is an astonishingly tight correlation between [X/Fe] and stellar age with amplitudes up to 0.20 dex over 
an age interval of eight Gyr in contrast to the lack of correlation between [Fe/H] and age. While [Mg/Fe] increases with age, the 
^-process element yttrium shows the opposite behavior meaning that [Y/Mg] can be used as a sensitive chronometer for Galactic 
evolution. [Na/Fe] and [Ni/Fe] are not well correlated with stellar age, but define a tight Ni-Na relation similar to that previously 
found for more metal-poor stars albeit with a smaller amplitude. Furthermore, the C/O ratio evolves very little with time, although 
[C/Fe] and [0/Fe] change by ~0.15 dex. 

Conclusions. The dependence of [X/Fe] on stellar age and the [Ni/Fe]- [Na/Fe] variations complicate the use of the [X/Fe]-rc relation 
as a possible signature for the existence of terrestrial planets around stars. The age trends for the various abundance ratios provide 
new constraints on supernovae yields and Galactic chemical evolution, and the slow evolution of C/O for solar metallicity stars is of 
interest for discussions of the composition of exoplanets. 

Key words. Stars: abundances - Stars: fundamental parameters - Stars: solar-type - (Stars:) planetary systems - Galaxy: disk - 
Galaxy: evolution 


1. Introduction 

In an important paper based on very precise determinations of 
abundances in the Sun and 11 solar twin stars, Melendez et al. 
(12009b found that the Sun has a lower refractory to volatile el¬ 
ement ratio than most of the solar twins. They suggest that this 
may be explained by the sequestration of refractory elements in 
terrestrial planets and that the slope of abundances relative to 
iron, [X/FejQ, as a function of elemental condensation tempera¬ 
ture, Tq, in a solar-composition gas (Lodders 12003b can be used 
as a signature of the existence of terrestrial planets around a star. 
This idea was supported by Ramirez et al. ( 12009b . who in a study 
of abundances of 64 solar-type stars find significant variations of 


* Based on data products from observations made with ESO 
Telescopes at the La Silla Paranal Observatory under programs given 
in Table [T] Tables 2 and 3 are available in electronic form at 
http://WWW.aanda.org 

* Eor two elements, X and Y, with number densities Nx and Ny, 
[X/Y] = log(Ax/AY)..ar - logfAx/Aylsun. 


the [X/Fej-Tc slope with only about 15% of the stars having a 
slope similar to that of the Sun. 

The works of Melendez et al. (12009b and Ramirez et al. 
(12009b have triggered several papers investigating relations be¬ 
tween [X/Fej-Tc slopes and existence of planets around stars 
(Gonzalez et al. 120101 Gonzalez Hernandez et al. 120101120131 
Schuler et al. 1201 11 Adibekyan et al. 120141 Maldonado et al. 
12015b . The results obtained are not so conclusive, partly because 
the stars included span larger ranges in T^ff, log g, and [Fe/H] 
than the solar twins in Melendez et al. (12009b making the abun¬ 
dance ratios less precise and introducing a dependence of the 
derived [X/Fej-Tc slopes on possible corrections for Galactic 
evolution effects. Thus, Gonzalez Hernandez et al. (12013b ob¬ 
tain both positive and negative slopes for ten stars with detected 
super-Earth planets, and Adibekyan et al. (12014b find evidence 
that the slope depends on stellar age. Furthermore, Schuler et 
al. (1201 lb show that the [X/Fej-Tc slope depends critically on 
whether one consider refractory elements only (Tq > 900 K) or 
include also volatile elements with Tq < 900 K. 
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Sequestration in terrestrial planets is not the only possible 
explanation of the low abundances of refractory elements in the 
Sun. Onehag et al. (120141 1 found that the [X/FeJ-Tc trend for 
solar-type stars belonging to the open cluster M 67 agrees with 
the trend for the Sun and suggest that the gas of the proto¬ 
cluster was depleted in refractory elements by formation and 
cleansing of dust before the stars formed. According to this, the 
Sun may have formed in a dense stellar environment contrary 
to most of the solar twins. Similar scenarios are discussed by 
Gaidos (I2015I I. who suggests that [X/Fel-Fc correlations can be 
explained by dust-gas segregation in circumstellar disks rather 
than formation of terrestrial planets. 

Additional information on the connection between abun¬ 
dance ratios in stars and the occurrence of planets may be ob¬ 
tained from studies of wide binaries for which the two compo¬ 
nents have similar and \ogg values. Dust-gas segregation 
before planet formation is likely to affect the chemical com¬ 
position of the components in the same way. In the case of 16 
CygA/B, for which the secondary component has a detected 
planet of at least 1.5 Jupiter mass (Mj), there is a difference 
[Fe/H](A-B)^ 0.04 dex (Ramirez et al. 1201II) and a difference 
in the [X/FeJ-Fc slope for refractory elements (Tucci Maia et al. 
12014b . which may be explained by accretion of elements in the 
rocky core of the giant planet around 16 CygB. Furthermore, 
the two components of XO-2, which are hosts of different plan¬ 
ets, have an abundance difference of ~ 0.05 dex for refractory 
elements (Teske et al. 120151 Damasso et al. 12015b . On the other 
hand, Liu et al. (12014b did not detect any abundance differences 
between the two components of HAT-P-1, for which the sec¬ 
ondary has a 0.5Mj planet. 

Recent precise determinations of abundances in four so¬ 
lar twin stars; HD101364 (Melendez et al. l20T2b . HIP 102152 
(HD197027) (Monroe et al. l20T3l l. 18Sco (HD 146233) 
(Melendez et al. l2014ab . and HD218544 (Melendez et al. l2014bb 
show that the hrst three mentioned stars have a similar depletion 
of refractory elements as the Sun, whereas 18 Sco has a [X/Fe]- 
Tq slope for elements from carbon to zinc similar to the majority 
of solar twin stars. Neutron capture elements in 18 Sco are, how¬ 
ever, signihcantly more abundant than lighter elements with the 
same condensation temperature, and in HIP 102152 several ele¬ 
ments (N, Na, Co, and Ni), have abundance ratios with respect 
to Fe that clearly fall below the [X/Fe]-Fc trend for the other el¬ 
ements. These anomalies suggest that the [X/Fej-Fc relation is 
not so well dehned as assumed in some of the works cited above, 
but that star-to-star variations in abundance ratios related to e.g. 
incomplete mixing of nucleosynthesis products or chemical evo¬ 
lution play a role. 

In order to make further studies of [X/Fe]-Fc relations, I 
have derived very precise abundance ratios (cr[X/Fe] ~0.01 dex) 
for a sample of 21 solar twin stars for which high-resolution 
HARPS spectra with signal-to-noise (S/N) ratios ranging from 
600 to more than 2000 are available. For the Sun a HARPS 
spectrum with S IN ^ 1200 of reflected sunlight from Vesta is 
used. The sample is not well-suited for studying relations be¬ 
tween abundance ratios and occurrence of planets, because only 
three of the stars are known to be hosts of planets. On the other 
hand, it is possible to study how well abundance ratios are cor¬ 
related with Fc and to see if [X/Fe] depends on other parameters 
such as stellar age. 

2. Observational data 

Based on the precise effective temperatures, surface gravities and 
metallicities derived by Sousa et al. (12008b from HARPS spec¬ 


tra of FGK stars, a sample of solar twins was selected to have 
parameters that agree with those of the Sun within +100K in 
Feff, +0.15 dex in log g, and ±0.10 dex in [Fe/H]. Listed in Table 
[1] 21 of these stars have HARPS spectra with S/N > 600 after 
combination of spectra available in the ESO Science Archive. 

The HARPS spectra have a resolution of R ^ 115 000 (Mayor 
et al. 12003b . Individual spectra for a given star were combined 
after correction for Doppler shifts and then normalized with the 
IRAF continuum task using cubic splines with a scale length of 
~ 100 A in the blue spectral region (3800 - 5300 A) and ~ 15 A in 
the red region (5350 - 6900 A), for which the archive spectra are 
more irregular. The S/N values given in Table [T] were estimated 
from the rms variation of the flux in continuum regions for the 
4700 - 6900 A spectral region, which contains the spectral lines 
applied in this study. By including spectra obtained from 2011 
to 2013 (programs 183.C-0972 and 188.C-0265), the S/N is in 
most cases signihcantly higher than that of the HARPS spectra 
applied in previous studies of solar twin stars (e.g., Gonzalez 
Hernandez et al. 12010112013b . For the stars in Table [1] the S/N 
ranges from 600 to 2400 with a median value of 1000. 

The IRAF splot task was used to measure equivalent widths 
(EWs) of spectral lines by Gaussian htting relative to pseudo¬ 
continuum regions lying within 3 A from the line measured. 
These regions do not necessarily represent the true continuum, 
but care was taken to use the same continuum windows in all 
stars, so that differences in EW between stars are precisely mea¬ 
sured given that the same instrument and resolution were applied 
for obtaining spectra. In this way differential abundances may 
be obtained with a precision better than 0.01 dex as shown by 
Bedell et al. (Iml . 

In order to achieve high-precision abundances relative to the 
Sun, it is very important that a solar flux spectrum is obtained 
with the same spectrograph as applied for the solar twin stars 
(Bedell et al. 12014b . As listed in Table [1] several such HARPS 
spectra are available based on reflected light from asteroids and 
the Jupiter moon Ganymede. As the prime reference spectrum 
for the Sun, the Vesta spectrum with S/N ^ 1200 is adopted, 
but it is noted that the Ceres and Ganymede spectra agree very 
well with that of Vesta. Equivalent widths measured in the Vesta 
spectrum and the combined spectrum of Ceres and Ganymede 
(S /N ^ 700) have an rms difference of 0.35 mA only. 

To illustrate the high quality of the spectra applied and the 
differential way of measuring EWs, Eig.[T] shows a comparison 
of the Vesta spectrum with that of HD 27063, which turns out 
to have a very similar effective temperature (Feff= 5779 K) and 
surface gravity (log g= 4.47) as the Sun. As seen from the differ¬ 
ence between the spectra in the lower panel of Eig. [1] lines cor¬ 
responding to refractory elements are somewhat stronger in the 
spectrum of HD 27063 than in the solar spectrum, whereas the 
line of the volatile element carbon at 5380.3 A is a bit stronger 
in the solar spectrum. HD 27063 turns out to be a star with a 
high refractory to volatile ratio compared to the Sun; i.e. [Ee/H]= 
0.064 and [C/H] = -0.018. 

Oxygen is an important volatile element in addition to car¬ 
bon. Unfortunately, the HARPS spectra do not cover the Oi 
triplet at 7774 A. Instead, the forbidden oxygen line at 6300.3A 
was used to derive [O/H]. Eor six stars telluric O 2 lines happen 
to disturb the [Oi] line in some of the spectra available. Such 
spectra were excluded before making a combination of spectra 
to be used in the [O i] region. This decreased the S/N of the [O i] 
spectra somewhat, but in all cases it is around or above 500. 
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Table 1. ESO/HARPS observing programs and S/N of combined 
spectra. 


Star 

ESO/HARPS program numbers 

S/N 

HD 2071 

072.C-0488, 

183.C-0972, 

188.C-0265 

1050 

HD 8406 

072.C-0488, 

183.C-0972 


700 

HD 20782 

072.C-0488, 

183.C-0972 


1250 

HD 27063 

072.C-0488 



850 

HD 28471 

072.C-0488, 

183.C-0972 


700 

HD 38277 

072.C-0488, 

183.C-0972 


900 

HD 45184 

072.C-0488, 

183.C-0972 


2400 

HD 45289 

072.C-0488, 

183.C-0972, 

188.C-0265 

1250 

HD 71334 

072.C-0488, 

183.C-0972, 

188.C-0265 

900 

HD 78429 

072.C-0488 



1200 

HD 88084 

072.C-0488, 

183.C-0972 


800 

HD 92719 

072.C-0488, 

183.C-0972 


1150 

HD 96116 

072.C-0488, 

183.C-0972, 

188.C-0265 

600 

HD 96423 

072.C-0488, 

183.C-0972, 

188.C-0265 

1200 

HD 134664 

072.C-0488 



850 

HD 146233 

072.C-0488 



850 

HD 183658 

072.C-0488, 

183.C-0972, 

188.C-0265 

1000 

HD 208704 

072.C-0488, 

183.C-0972, 

188.C-0265 

1100 

HD 210918 

072.C-0488 



1300 

HD 220507 

072.C-0488, 

183.C-0972, 

188.C-0265 

1100 

HD 222582 

072.C-0488, 

183.C-0972 


1000 

Vesta 

088.C-0323 



1200 

Ceres 

082.C-0357, 

060.A-9036, 

060.A-9700 

500 

Ganymede 

060.A-9036, 

060.A-9700 


450 



Fig. 1. The solar flux HARPS spectrum around the C i line at 
5380 A obtained via reflected light from Vesta (thin, red line) in 
comparison with the HARPS spectrum of HD 27063 (thick, blue 
line). The lower panel shows the difference (HD 27063 - Sun) 
between the two spectra. 


3. Model atmosphere analysis 

The Uppsala EQWIDTH program, which is based on the as¬ 
sumption of local thermodynamic equilibrium (LTE), was used 
to calculate equivalent widths as a function of element abun¬ 
dance for plane parallel (ID) MARCS model atmospheres 
(Gustafsson et al. 120081 1 representing the stars. Interpolation to 
the observed EWs then yields the LTE abundances for the vari¬ 
ous elements. 


3.1. Spectral lines 

The spectral lines used are listed in Table |2] As seen, they fall 
in the spectral region 4700 - 6900 A, where the S/N of HARPS 
spectra is higher and line blending less severe than in the blue 
region. Strong lines, such as the NaD lines and the Mgb triplet, 
and lines significantly blended by stellar or telluric lines were 
avoided except for the [O i] 46300.3 line, which is blended by a 
Ni I line. Eor C, Na, Mg, and Al two lines are available for each 
element. Heavier elements have at least three lines and for Ti, Cr, 
and Ee both neutral and ionized lines are available. In the case 
of Ee, only lines having EWsun < 70 mA were included. This 
results in a list of 47 Eei lines with excitation potentials, ;t'exc, 
ranging from 0.9 to 5.1 eV and nine Eeii lines, which are used 
to determine effective temperatures and surface gravities of the 
stars as described in Sect. 13.21 

As the analysis is made differentially to the Sun line by line, 
the gf/-values of the spectral lines cancel out. Doppler broaden¬ 
ing due to microturbulence is specified by a depth-independent 
parameter, ^turb- Collisional broadening of Cat, Tii, Cri, Eei, 
Ee It, and Ni i lines caused by neutral hydrogen and helium atoms 
is based on quantum mechanical calculations (Anstee & O’Mara 
IT9951 Barklem & 0’ Mara[T998l Barklem et al. lIOOOl Barklem & 
Asr)elund-Johansson l2005t . Eor the remaining lines, the Unsold 
(119551) approximation with an enhancement factor of two was 
applied. If instead an enhancement factor of one is adopted, the 
differential abundances are not changed significantly due to the 
fact that the damping wings of the lines only contribute a small 
fraction of the equivalent width. 

The derivation of oxygen abundances impose a particular 
problem, because the [Oi] 46300.3 line is blended by a Ni i line 
with nearly the same wavelength (Allende-Prieto et al. 120011 ). 
Based on the MARCS solar model, loggf - -2.11 (Johansson 
et al. 120031 ). and a solar nickel abundance of A(Ni)Q = 6.15 (as 
determined from 17 weak Nit lines listed in Scott et al. 120091 ). 
the equivalent width of the Ni i line is calculated to be 1.7 mA. 
After subtracting this value from the measured EW of the [O i] 
- Nit blend (5.3mA), one gets EW([Oi])o = 3.6mA. Eor the 
stars, the measured EWs of the [Oi] - Ni i blend were corrected 
in the same way by subtracting the EW of the Ni i line as calcu¬ 
lated for the model atmosphere of the star and the [Ni/H] value 
derived. 

3.2. Stellar parameters 

Adopting solar parameters, refr = 5777K, log 0 = 4.438, and 
^tuib = 1.0 km s“', the corresponding parameters of the stars were 
determined by requesting that [Ee/H] has no systematic depen¬ 
dence on Xexc and EW of the lines and that the mean [Ee/H] val¬ 
ues derived from Ee i and Ee ii lines, respectively, are equal. The 
slope of [Ee/H] as a function of excitation potential is sensitive 
to the effective temperature and the slope of [Ee/H] versus EW 
depends on the microturbulence. The difference of [Ee/H] de¬ 
rived from Ee i and Ee it lines depends on the surface gravity via 
its effect on the electron pressure in the stellar atmospheres, but 
it is also affected by the iron and a-element abundances, because 
these elements are important electron donors. Therefore, it it is 
necessary to make a number of iterations using MARCS models 
with different parameters Tefr, log g, [Ee/H], and [a/Ee] to ensure 
that the excitation and ionization balance of Ee lines are fulfilled. 
These models are obtained from two sets of MARCS models; the 
standard set with [a/Ee] = 0 at [Ee/H] > 0 and [a/Ee] rising lin¬ 
early to [a/Ee]= 0.4 when [Ee/H] decreases from 0.0 to -1.0, 
and the a-enhanced set with [a/Ee] = 0.4 at all metallicities. 
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Exo. pol. [eV] 



4.6 


4.5 


x: 

O) 

g? 4.3 


4.2 


4.1 

4.1 4.2 4.3 4.4 4.5 4.6 

logs (spec.) 

Fig. 3. Photometric gravities based on Hipparcos parallaxes as a 
function of spectroscopic gravities determined from the Fe ion¬ 
ization balance. 



Fig. 2. [Fe/H] derived from lines in the spectrum of HD 183658 
as a function of excitation potential (upper panel) and “normal¬ 
ized” equivalent width: EWamm. - EW ■ 6000.0/A [A] (lower 
panel). Open circles refer to Fei lines and filled (red) circles 
to Fe II lines. The full drawn lines show linear fits to the Fe i 
data for the derived parameters of the star. The dashed lines 
show changes in the slopes if Tefr is decreased by 30 K (upper 
panel) and the microturbulence is decreased by 0.1 km s * (lower 
panel). 

Hence, it it possible to interpolate to a MARCS model with the 
[a/Fe] valu^ of a star. This is of some importance for three a- 
enhanced solar twins in this paper with [a/Fe] ^ 0.08. If models 
with [a/Fe] = 0.0 were adopted for these stars, the derived ef¬ 
fective temperature and gravity would change by ATeff ^ -8 K 
and Alog g ^ -1-0.025 dex, whereas the effect on [Fe/H] is only 
-0.003 dex. 

As a typical example. Fig. |2] shows [Fe/H] as a function of 
Xexc and EW obtained with the final model of HD 183658 (Teff, 
log g, [Fe/H], [a/Fe]) = (5809 K, 4.402, h- 0.035, 0.007). The rms 
scatter of [Fe/H] is 0.007 dex for Fe i lines and 0.011 dex for Fe ii 
lines. Based on the 1-sigma errors of the slopes of [Fe/H] ver¬ 
sus ;^fexc and EW and the errors of the average Ee abundances 
derived from Ee i and Ee ii lines, respectively, the error analy¬ 
sis method described by Epstein et al. (120101 1 and Bensby et al. 
(12014b has been used to derive the errors of Teir, logg, [Ee/H], 
and ^turb- The estimated errors are nearly the same for all stars, 
i.e., cr(refr) = 6K, (T(\ogg) - 0.012 dex, cr([Ee/H]) = 0.006 dex, 
and cr(^tm-b) = 0.017km s“*. 

In the derivation of the model atmosphere parameters, small 
differential non-LTE corrections for the Fe i lines according to 
Lind et al. (12012b have been included. The largest non-LTE ef¬ 
fects occur at the lowest gravities of the solar twin sample, i. e., 
log g ~ 4.25; [Ee/H] is increased by 0.003 to 0.004 dex and log g 
by 0.007 to 0.009 dex relative to the LTE values. Changes in 
are negligible, i. e., ranging from -i-2 to -3 K. 

The small errors quoted are only obtainable for samples of 
solar twin stars that can be analyzed differentially with respect 

2 Defined as [a/Fe] = \ ■ ([Mg/Fe] -b [Si/Fe] + [Ca/Fe] + [Ti/Fe]) 
in this paper. 


to the Sun using the same list of Fe lines. Still, one may wonder 
if these errors are really realistic. To test this, I have compared 
with parameters recently determined for 88 solar twin stars by 
Ramirez et al. (I2014bb based on spectra obtained with the 6.5 m 
Clay Magellan Telescope and the MIKE spectrograph. The spec¬ 
tra have S IN > 400 and resolutions of R = 83 000 (65 000) in 
the blue (red) spectral regions. A set of 91 Fe i and 19 Fe ii lines 
were analyzed with MARCS models to derive atmospheric pa¬ 
rameters with errors similar to those quoted above. For 14 stars 
in common with the present paper I get the following average 
differences (Ramirez - this paper) and rms deviations: ATeft = 
0 + lOK, Alog 0 = 0.002 + 0.020, and A[Fe/H] = 0.000 + 0.014. 
The rms deviations for Teir and log g are fully explainable by the 
errors quoted in the two works, whereas the rms deviation for 
[Fe/H] is a bit larger than expected. 

The parameters derived in this paper are also in good agree¬ 
ment with those derived from HARPS spectra by Sousa et al. 
(12008b based on an extensive list of 236 Fe i and 36 Fe ii lines 
analyzed by using a grid of Kurucz model atmospheres (Kurucz 
11993b . The average differences (Sousa - this paper) and rms de¬ 
viations for the 21 solar twins in the present paper are: ATeit = 
-1+8K, Alog 0 = 0.018+0.033,andA[Fe/H] = -0.003 + 0.009. 
The deviations in and [Fe/H] are compatible with the small 
errors quoted above, but the rms deviation in log g is larger than 
expected; stars with \ogg < 4.35 turn out to have systematically 
higher log g values in Sousa et al. (12008b than in the present pa¬ 
per. These deviations are not seen when comparing with logg 
values from Ramirez et al. (I2014bb . 

The spectroscopic surface gravities have also been compared 
to photometric gravities determined from the relation 

n ]\/f TpfF 

log = log — + 4 log --+ 0.4(Mboi - Mboi, 0 ), (1) 

go Mq Teff,© 

where M is the mass of the star and Mboi the absolute bolomet- 
ric magnitude. Hipparcos parallaxes (van Leeuwen 12007b were 
used to derive absolute visual magnitudes from V magnitudes 
based on Stromgren photometry (Olsen 11983b . and bolometric 
corrections were adopted from Casagrande et al. (12010b . Stellar 
masses were obtained by interpolating in the luminosity - log Teir 
diagram between the Yonsei - Yale evolutionary tracks of Yi et al. 
(12003b : see Nissen & Schuster (12012b for details. The estimated 
error of log g (phot.) ranges from 0.02 to 0.04 dex depending pri- 
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marily on the relative error of the parallax. A seen from Fig. [3] 
there is a satisfactory agreement between the two sets of gravi¬ 
ties. The weighted mean deviation (log o'(phot.) - log o'(spec.)) 
is 0.005 dex, and the rms deviation is as expected from the error 
bars except for one star, HD 38277, which has a 3-sigma devi¬ 
ation. Its lower log 0 (phot.) could be explained if the star is a 
spectroscopic binary, but there is no indication of a component 
in the HARPS spectra. 

3.3. Non-LTE and 3D effects 

Although the solar twin stars have atmospheric parameters close 
to those of the Sun, effects of deviations from LTE may be im¬ 
portant considering the very high precision we are aiming at. 

The detailed calculations by Lind et al. (I2012I I show that Fe ii 
lines are not affected by departures from LTE, but as mentioned 
in Sect. 13.21 there is a small non-LTE effect on Fe abundances 
derived from Fe i lines. Using the IDE program made available 
by Lind et al. to calculate non-LTE corrections for the actual 
set of Fe i lines applied, it is found that the largest differential 
correction, A[Fe/H] = 0.004, occurs for stars having gravities 
about 0.15 dex below the solar gravity. Calculations of non-LTE 
corrections as a function of Teff, log g, and [Fe/H] for the lines 
of CI (Takeda & Hon dal2005l l. Mg i (Zhao & Gehren l2000t . S i 
and Zn i (Takeda et al. l2005t . and Ca i (Mashonkina et al. 12007b . 
suggest similar small differential corrections as in the case of the 
Fei lines; hence the non-LTE effect on the corresponding abun¬ 
dance ratios, [X/Fe], is negligible, and the LTE ratio is therefore 
adopted. For the two Nai lines applied, 46154.2 and 46160.7, 
the calculations of Lind et al. (1201 lb . indicate similar numerical 
corrections as in the case of Fe i lines, but with opposite sign. 
Therefore, non-LTE corrections have been applied when deter¬ 
mining [Na/Fe] ratios; the most extreme correction is -0.008 dex 
for HD 38277. 

For the lines of Al i and Si i there seems to be no detailed 
calculations of non-LTE effects as a function of Teff, log g, and 
[Fe/H] for solar type stars, but the corrections for the Sun itself 
are on the order of 0.01 dex only (Baumiiller & Gehren lT996l Shi 
et al. 120081 Scott et al. I2015bb . This suggests that differential 
corrections for the Al and Si abundances can be neglected. It 
is probably also the case for Ni (Scott et al. I2015ab . although 
extensive non-LTE calculations for Ni i lines are still missing. 

The statistical equilibrium calculations by Bergemann 
(1201 lb and Bergemann & Cescutti (12010b show that the Ti ii and 
Cm lines applied in this paper are formed in LTE in the solar 
atmosphere, whereas there are significant non-LTE effects for 
the Til and Cri lines. We may investigate this empirically by 
comparing abundances derived from neutral and ionized lines, 
respectively. As seen from Fig. |4] there appears to be a trend of 
[Ti i/Ti ii] as a function of log g suggesting that solar twin stars 
with the lowest and highest gravities have differential non-LTE 
corrections on the order of +0.01 dex for Ti abundances derived 
from Ti i lines. In view of this, the Ti abundances derived from 
Tin lines will be adopted although only three such lines are 
available compared to 11 Tii lines. For Cr there is no signifi¬ 
cant trend of [Cri/Crii] versus loggr and the rms scatter is only 
0.008 dex. In this case, the Cr abundance derived from Cr i lines 
is adopted. 

The oxygen abundances derived from the forbidden [Oi] 
line at 6300 A are not affected by departures from LTE (e.g., 
Kiselman 11993b . This line is , however, sensitive to 3D hydro- 
dynamical (granulation) effects on the temperature structure of 
the stellar atmosphere (Asplund 12005b . According to Nissen et 



logs 

Fig. 4. The ratio of Ti abundances derived from Ti i and Ti ii 
lines, respectively, as a function of log g (upper panel) and the 
corresponding ratio for Cr (lower panel). 

al. (12002b . this leads to 3D corrections that depend on effective 
temperature and metallicity. For the present sample of solar twin 
stars the largest differential effect on [O/H] is 0.007 dex. For the 
other elements, one would expect 3D effects to be less signifi¬ 
cant, because the lines applied are formed deeper in the atmo¬ 
spheres, but this need to be confirmed by 3D model calculations 
preferentially with non-LTE effects included. 

The impact of magnetic fields on the temperature structure 
of stellar atmospheres is another potential problem. Magneto¬ 
hydrodynamic simulations of the solar atmosphere (Fabbian et 
al. 12010112012b show effects on the strength of Fe i lines induced 
primarily by “magnetic heating” of the upper layers. For an av¬ 
erage vertical magnetic field strength of 100 G (Trujillo Bueno 
et al. 12004b . the effect on the derived Fe abundance ranges from 
0.01 dex for weak lines formed in deep layers to 0.08 dex for 
stronger lines formed in the upper part of the photosphere. These 
results have recently been confirmed by Moore et al. (12015b . who 
in addition find a correction on the order of + 0.02 dex for the 
oxygen abundance derived from the [O i] 46300 line if a vertical 
mean field of 80 G is imposed. A similar correction is suggested 
by the recent MHD simulations of Fabbian & Moreno-Insertis 
(l20T5Tl . 

These simulations raise the question if star-to-star differ¬ 
ences of the magnetic field strength lead to variations of abun¬ 
dance ratios among solar twin stars. If so, the ratio of Fe abun¬ 
dances determined from Fe i and Fe ii lines (used for determining 
the gravity parameter) would also be affected, because of differ¬ 
ent depth of formation for neutral and ionized lines. As discussed 
in Sect. 13.21 there is, however, good agreement between spectro¬ 
scopic gravities and gravities determined via Hipparcos paral¬ 
laxes, which suggests that differences in magnetic field strength 
do not play a major role in the context of abundance variations 
among solar twin stars. Still, this is a problem that deserves fur¬ 
ther studies. 

3.4. Abundance ratios and uncertainties 

The derived abundance ratios, [X/Fe], are given in Table [3] to¬ 
gether with the atmospheric parameters of the stars. Due to the 
systematic deviation between log g values derived in this paper 
and those of Sousa et al. (12008b . four stars have gravities below 
the selection limit of log g - 4.29. 
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Table 4. Uncertainties of abundance ratios. 



^lines 

^atm.par. 

[dex] 

C EW 

[dex] 

^^adopted 

[dex] 

[C/Fe] 

2 

0.009 

0.009 

0.013 

[0/Fe] 

1 

0.008 

indv. 

indv. 

[Na/Fe] 

2 

0.002 

0.004 

0.007 

[Mg/Fe] 

2 

0.002 

0.010 

0.011 

[Al/Fe] 

2 

0.002 

0.006 

0.008 

[Si/Fe] 

10 

0.004 

0.002 

0.007 

[S/Fe] 

4 

0.009 

0.012 

0.016 

[Ca/Fe] 

7 

0.002 

0.003 

0.006 

[Ti/Fe] 

3 

0.007 

0.007 

0.010 

[Cr/Fe] 

10 

0.001 

0.003 

0.006 

[Ni/Fe] 

14 

0.001 

0.003 

0.006 

[Zn/Fe] 

3 

0.004 

0.009 

0.011 

[Y/Fe] 

3 

0.007 

0.005 

0.010 


The estimated 1-sigma errors of the differential abundance 
ratios are given in Table |4] Column three lists the error cor¬ 
responding to the uncertainty of the atmospheric parameters 
(o-CTeff) = +6K, cr(log0) = +0.012, cr([Fe/H]) = 0.006 dex 
and cr(^turb) = 0.017kms“'). For some of the abundance ratios, 
this error is very small, because the derived abundances of ele¬ 
ment X and Fe have almost the same dependence on the param¬ 
eters. Column four lists the error arising from the measurement 
of equivalent widths, calculated as the line-to-line scatter of the 
abundances derived divided by VMines - 1, where Mines is the 
number of spectral lines applied. In the case of oxygen for which 
only the weak [Oi] /16300.3 line was used, c+ew was calculated 
by assuming that the error of the equivalent width measurement 
is cr{EW) - ±0.2 mA. Finally, I have taken into account that dif¬ 
ferential 3D effects may introduce errors in the abundance ratios 
up to 0.005 dex. Adding these error sources in quadrature leads 
to the total adopted errors listed in the last column of Tabled 


4. Stellar ages 

As discussed in Sect. |5] it is interesting to have information on 
stellar ages. Although the stars belong to the main sequence, rel¬ 
ative ages may be derived from the precise values of Teff and 
logg as shown in Fig. |5] where the stars are compared to a 
Yonsei-Yale set of isochrones (Yi et al. 120011 Kim et al. 120021) . 
Using the IDL program made available by Kim et al. to calculate 
isochrones for the actual combination of [Fe/H] and [a/Fe] for 
a given star, it is possible to interpolate to the stellar age and its 
internal 1-sigma error. The results are given in Table |5] together 
with stellar masses derived from and log g. The internal error 
of the mass is on the order of iO.OlMsun- 

As seen from Fig.|6l the stellar ages agree very well with the 
values of Ramirez et al. (I2014bl l. who also determined ages from 
Yonsei-Yale isochrones in the Teff- log g diagram. The small sys¬ 
tematic deviation for the three oldest stars may be due to the fact 
that they have enhanced [(r/Fe] values ([a/Fe] ^ 0.08), which 
were taken into account when deriving their ages in the present 
paper, but apparently not by Ramirez et al. (I2014bl) . 



5900 5850 5800 5750 5700 5650 

n,, [K] 


Fig. 5. The Teff- log g diagram of solar twin stars in comparison 
with Yonsei-Yale isochrones corresponding to ages given to the 
right in Gyr. The isochrones refer to [Fe/H]= 0.0 and [a/Fe] = 
0.0, but 5 Gyr isochrones corresponding to [Fe/H] = -0.05 
(dashed line) and [a/Fe] = +0.1 (dotted line) are also shown. 
The position of the Sun with an age of 4.55 Gyr is marked. 



Fig. 6. Comparison of stellar ages derived in this paper with ages 
determined by Ramirez et al. (I2014bll . 

5. Discussion 

5.1. [X/Fe] - stellar age correlations 

As seen from Fig. |7] the dispersion of [X/Fe] at a given [Fe/H] 
is larger than expected from the estimated error bars except in 
the case of [Cr/Fe]. Three stars (HD 45289, HD 210918, and 
HD 220507) marked by open (red) squares stand out by having 
high [X/Fe] ratios of the alpha-capture elements O, Mg, Si, S, 
and Ti, and are also high in C, Al, and Zn. These stars proba¬ 
bly belong to the class of metal-rich a-enhanced stars first iden¬ 
tified by Adibekyan et al. (1201 11 1201211 as a population hav¬ 
ing [a/Fe] ^ +0.1 in contrast to normal thin-disk stars with 
[a/Fe] ^ +0.0. Bensby et al. (120141) confirm the existence of 
such metal-rich a-enhanced stars, and Haywood et al. (120131) 
show that they belong to a thick-disk sequence in the [a/Fe] 
- [Fe/H] diagram for which the ages range from ~ 13 Gyr at 
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Fig. 7. Abundance ratios [X/Fe] as a function of [Fe/H]. The error bars shown correspond to the errors listed in Table |4] Open (red) 
squares: a-enhanced stars. Filled circles: thin-disk stars 


[Fe/H] = -1 to ~ 8Gyr at [Fe/H] = +0.2. The ages of the 
three cr-enhanced stars found in this paper (8-9 Gyr) agree with 
this interpretation. Furthermore, two of the a-enhanced stars 
(HD45289 and HD 210918), have thick-disk kinematics accord¬ 
ing to the Galactic velocity components given in Adibekyan et 
al. (120121 1. The third one has typical thin-disk kinematics like the 
rest of the sample. 

Even if we exclude the three a-enhanced stars, there are still 
variations in the abundance ratios that cannot be explained by 
the estimated errors of the determinations. Interestingly, there 
is a clear correlation between [X/Fe] and stellar age for most 
of the elements as seen from Fig. |8] Linear fits ([X/Fe] - a + 
b ■ Age) have been obtained by a maximum likelihood program 
that includes errors in both coordinates. Coefficients a and b are 
given in Table |6] together with standard deviations and reduced 
chi-squares of the hts. As seen, [Mg/Fe], [S/Fe], [Ti/Fe], and 
[Cr/Fe] are very well correlated with age ~ 1). but [C/Fe], 
[O/Fe], [Al/Fe], [Si/Fe], [Cr/Fe], and [Zn/Fe] also show a clear 
dependence on age. Exceptions are [Na/Ee] and [Ni/Ee], which 
will be discussed below. 

Eigure |3 (lower panel) shows how astonishingly well 
[Mg/Ee] for thin-disk stars is correlated with age in contrast to 
the lack of correlation between [Ee/H] and age (upper panel). 
Given that Type la Supernovae (SNe) provide about two-thirds 
of Ee in the Galactic disk but very little Mg (Kobayashi et al. 


I 2 OO 6 I 1 . the age correlation of [Mg/Ee] can be explained by an in¬ 
creasing number of Type la SNe relative to the number of Type 
II SNe as time goes on. The small scatter in the [Mg/Ee]- age 
relation then indicates that nucleosynthesis products of Type la 
and Type II SNe are well mixed locally in the Galactic disk 
before new stars are formed. The large scatter in metallicity 
(cr[Ee/H] ^ 0.07 dex) at a given age may be explained by infall 
of metal-poor gas onto the Galactic disk followed by star for¬ 
mation before mixing evens out the chemical inhomogeneities. 
The infalling gas will mainly contribute hydrogen causing a de¬ 
crease of [Ee/H] but leave [Mg/Ee] unaffected. This scenario is 
discussed in detail by Edvardsson et al. (119931 1. who were hrst 
to discover the large spread in the age-metallicity relation of 
thin-disk stars in contrast to a small spread in [a/Eek but their 
precision of abundance ratios (typically ±0.05dex)[3 was not 
high enough to see a correlation between [a/Ee] and stellar age. 
Alternatively, the dispersion in [Ee/H] at a given age may be re¬ 
lated to orbital mixing of stars born at different galactocentric 
distances (Rg) provided that there is a strong radial gradient in 
[Ee/H] around the solar distance (Haywood et al. l2013ll . but then 
the very small dispersion in [X/Ee] at a given age requires that 
the star formation history is independent of Rg (Haywood et al. 

IMBll . 


^ The Edvardsson et al. sample is not limited to solar twin stars and 
abundances were derived from high-resolution spectra with S /N ~ 200. 
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Fig. 8. Abundance ratios [X/Fe] as a function of stellar age for the thin-disk stars. The lines show maximum likelihood linear hts to 
the data, with zero points and slope coefficients as given in Table |6l 


As seen from Fig. |8] and Table |6l the slope of [X/Fe] versus 
age differs from element to element. To some extent this may 
be explained in terms of differences in the relative contributions 
of Type la and Type II SNe to the various elements, but it can¬ 
not explain that the age-slopes for C, O, and A1 are significantly 
larger than that of Mg, because for all four elements, the yield 
of Type la SNe is negligible compared to the Type II SNe yield 
(e.g., Kobayashi et al. l2006l l. Hence, it seems that one has to in¬ 
clude an evolving initial mass function to explain the age trends 
seen in Fig. [8] Asymptotic giant branch (AGB) stars (Karakas 
1201 01 1 should also be considered, but cannot solve the problem; 
they contribute C and A1 (as well as Na) but no Fe, and hence 
increase [X/Fe] for these elements as a function of increasing 
time, i.e., decreasing stellar age. 

The near-constancy of [Ca/Fe] as a function of stellar age 
(see Fig. [8]l is another puzzling problem. According to yields of 
supernovae (Kobayashi et al. I2006t . one would expect [Ca/Fe] 
to increase with age in the same way as [Si/Fe] and [S/Fe]. A 
similar problem is encountered in studies of abundance ratios in 
early-type galaxies (Conroy et al. 1201411 : [Ca/Fe] is constant as a 
function of stellar velocity dispersion in contrast to a rising trend 
for the other a-capture elements. Mulchaey et al. (120141) suggest 
that this special behavior of Ca as well as the high Ca abundance 
in diffuse gas of clusters of galaxies is due to a new class of 
low luminosity supernovae exemplified by SN 2005E, for which 
the amount of synthesized Ca is 5-10 times greater than that of 


classical Type la SNe (Perets et al. 12010b . Perhaps these Ca-rich 
SNe are also important for the chemical evolution of the Galactic 
disk. 

Diffusion of elements (gravitational settling and radiative 
levitation) should also be considered when discussing abun¬ 
dances of stars as a function of time. Helioseismology and mod¬ 
els of the Sun indicate that the surface abundances of heavy el¬ 
ements (and helium) have decreased by ~ 10% relative to their 
initial values (Christensen Dalsgaard et al. 11996b . According to 
Turcotte & Wimmer-Schweingruber (12002b . the effect on [X/Fe] 
is, however, much smaller, i.e., < 0.004 dex for the elements dis¬ 
cussed in this paper. This means that diffusion can be neglected 
when discussing the age trends of [X/Fe]. 

The i-process element yttrium is not included in Fig. [8j but 
in contrast to most of the other elements, [Y/Fe] decreases with 
stellar age as previously found for barium (Edvardsson et al. 
119931 Bensby et al. 12007b . This is probably due to an increasing 
contribution of Y and Ba from low-mass (1-4 Msun) AGB stars 
as time go on (Travaglio et al. 12004b . As the slope of [Y/Ee] ver¬ 
sus age is quite steep {b - -0.033 dex Gyr^’ with opposite sign 
to that of [Mg/Ee] {b = -1-0.009 dex Gyr^'), [Y/Mg] becomes 
a sensitive indicator of age as seen from Eig. [10] A maximum 
likelihood fit to the data yields the relation 

[Y/Mg] = 0.175 (+0.011) - 0.0404 (+0.0019) Age [Gyr] (2) 

with;rL - 0.71. 
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Table 5. Stellar ages and masses. 


Star 

Age 

[Gyr] 

Mass 

[Mo] 

Note" 

HD 2071 

3.5 ±0.8 

0.97 


HD 8406 

4.2 ±0.8 

0.96 


HD 20782 

7.5 ±0.4 

0.97 

1.8Mj* 

HD 27063 

2.6 ±0.6 

1.04 


HD 28471 

7.0 ±0.4 

0.97 


HD 38277 

7.3 ±0.4 

1.01 


HD 45184 

2.7 ±0.5 

1.06 

0.04Mj‘’ 

HD 45289* 

8.5 ±0.4 

1.00 


HD 71334 

8.1 ±0.4 

0.94 


HD 78429 

7.5 ±0.4 

1.04 


HD 88084 

6.0 ±0.6 

0.96 


HD 92719 

2.7 ±0.6 

0.99 


HD 96116 

0.7 ±0.7 

1.05 


HD 96423 

6.0 ±0.4 

1.03 


HD 134664 

2.3 ±0.5 

1.07 


HD 146233 

3.8 ±0.5 

1.04 


HD 183658 

5.0 ±0.5 

1.03 


HD 208704 

6.9 ±0.4 

0.98 


HD 210918* 

8.5 ±0.4 

0.96 


HD 220507* 

8.8 ±0.4 

1.01 


HD 222582 

6.5 ±0.4 

1.01 

i.m/ 


Notes. a-enhanced star. Mass of detected planet in units of 
Jupiter’s mass. ® Jones et al. 120061 . Mayor et al. 120111 . Butler 
et al. 120061 . 


Table 6. Linear fits of [X/Fe] as a function of stellar age. 



a 

[dex] 

b 

10“^ dex Gyr-‘ 

o-[X/Fe] “ 
[dex] 

yL 

[C/Fe] 

-0.110 

+ 

0.011 

±13.9 

±2.0 

0.025 

2.7 

[O/Fe] 

-0.098 

+ 

0.015 

±16.1 

±2.8 

0.030 

1.9 

[Na/Fe] 

-0.090 

+ 

0.035 

±9.8 

±6.3 

0.039 

20.6 

[Mg/Fe] 

-0.037 

+ 

0.005 

±8.9 

±0.9 

0.011 

0.8 

[Al/Fe] 

-0.074 

+ 

0.011 

±16.7 

± 1.9 

0.020 

3.6 

[Si/Fe] 

-0.028 

+ 

0.005 

±5.7 

±0.9 

0.012 

2.4 

[S/Fe] 

-0.058 

+ 

0.007 

±7.0 

± 1.3 

0.015 

0.8 

[Ca/Fe] 

±0.020 

+ 

0.004 

-0.6 

±0.8 

0.010 

2.7 

[Ti/Fe] 

-0.013 

+ 

0.005 

±6.3 

±0.8 

0.010 

0.9 

[Cr/Fe] 

±0.014 

+ 

0.003 

-2.6 

±0.6 

0.007 

1.4 

[Ni/Fe] 

-0.044 

+ 

0.011 

±4.8 

±2.0 

0.021 

10.7 

[Zn/Fe] 

-0.079 

+ 

0.008 

±11.7 

± 1.4 

0.018 

2.0 

[Y/Fe] 

±0.146 

+ 

0.011 

-33.0 

±2.0 

0.023 

1.6 

Notes. Standard deviation of [X/Fe] for the linear fit. 


For elements having a condensation temperatures markedly 
different from that of Fe, the variations in the [X/FeJ-Tc slope 
discussed in Sect. l5.2l contribute to the scatter in the [X/Fe]-age 
relation. In particular, this effect may explain why the scatter 
in [C/Fe] and [O/Fe] is larger than expected from the estimated 
errors of the determinations, = 2.7 and 1.9, respectively. 

[C/O], on the other hand, is not much affected by the [X/Fe]- 
Tq slope variations, because the two elements have similar low 
condensation temperatures. 



0 2 4 6 8 10 

Age [Gyr] 


Fig. 9. [Fe/H] and [Mg/Fe] versus stellar age. Thin-disk stars are 
shown with filled circles and a-enhanced stars with open (red) 
squares. The line in the lower panel corresponds to the fit given 
in Table |6l 



Age [Gyr] 

Fig. 10. [Y/Mg] versus stellar age for thin-disk stars. 

Recently, the C/O ratioQ in solar-type stars has been much 
discussed, mainly because its value may influence the compo¬ 
sition of planets formed in circumstellar disks (e.g.. Bond et al. 
120101 Fortney 120121 Gaidos 120151 . In some studies (Delgado 
Mena et al. 120101 Petigura & Marcy 1201 II I. C/O was found to 
range from ~ 0.4 to i 1.0, i.e., up to a factor of two higher 
than the solar ratio, C/O© = 0.55 + 0.10 (Asplund et al. 120091 


C/O is defined as Nq/Nq, where Nq and No are the number densities 
of carbon and oxygen nuclei, respectively, and should not be confused 
with the solar-normalized logarithmic ratio, [C/O]. 
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Fig. 11. The C/O ratio as a function of stellar age. The line shows 
the evolution of C/O (normalized to a solar ratio of 0.55) accord¬ 
ing to the chemical evolution calculations of Gaidos (12015b . 

Caffau et al. 1201 lb . This has led to speculations about the exis¬ 
tence of exoplanets consisting of carbides and graphite instead 
of Earth-like silicates (Bond et al. 12010b . Alternative determi¬ 
nations of C/O in solar-type stars (Nissen 120131 Nissen et al. 
120141 Teske et al. 12014b show, however, only a small scatter in 
C/O and an increasing trend as a function of [Fe/H] reaching 
C/O ^ 0.8 at [Fe/H] 0.4. The scatter of C/O for the solar 
twin stars is also small, i.e., +0.035 only around a mean value of 
<C/0> = 0.52 (see Fig.fTTb. Furthermore, there is not much evo¬ 
lution of C/O as a function of time despite the fact that both C/Fe 
and O/Fe evolve significantly. This agrees very well with the re¬ 
cent Galactic chemical evolution calculations of Gaidos (12015b 
as seen from Fig. [TT] In his model, the dominant contribution of 
C and O comes from massive stars, and C/O has only a moderate 
rise during the first ~ 5 Gyr due to an increasing contribution of 
carbon from stars of lower mass. 

As mentioned above, [Na/Fe] and [Ni/Fe] are not well cor¬ 
related with stellar age; according to Tablej^the linear fits have 
Tred ^ Curiously, there is, however, a tight correlation be¬ 
tween [Na/Fe] and [Ni/Fe] as shown in Fig. [12] A maximum 
likelihood linear fit with errors in both coordinates provides the 
relation 

[Ni/Fe] = 0.002 (±0.002) + 0.580 (±0.032) [Na/Fe] (3) 

with a standard deviation cr[Ni/Fe] = 0.008, which is close to 
that expected from the estimated errors of [Na/Fe] and [Ni/Fe], 
i. e., ±0.007 and ±0.006 dex, respectively. 

Apparently, correlated variations of [Ni/Fe] and [Na/Fe] 
among solar-metallicity stars have not been noted before. This is 
understandable given that the amplitude is small, i.e., ~0.07dex 
in [Ni/Fe]; it requires very high abundance precision to see 
the variations^. A Ni-Na relation has, however, been found for 
thick-disk and halo stars with -1.6 < [Fe/H] - 0.4 (Nissen & 
Schuster 1199711201011201 lb and for stars in dwarf spheroidal 
galaxies in the same metallicity interval (Venn et al. 120041 
Fetarte l70101 Femasle [2014b . For these metal-poor stars, the am¬ 
plitude of the variations is larger (~ 0.3 dex in [Ni/Fe]), but the 
slope of the [Ni/Fe] - [Na/Fe] relation is about the same as found 
for the solar-metallicity stars, i. e. A[Ni/Fe]/A[Na/Fe] ^ 0.5. 
Furthermore, the lowest Ni/Fe and Na/Fe ratios are found among 

^ For the solar twin stars studied, abundance ratios derived from 
HARPS spectra by Neves et al. ( 12009b show in fact correlated varia¬ 
tions in [Ni/Fe] and [Na/Fe] similar to Fig.[T2| but it was not noted in 
that paper nor in subsequent papers based on HARPS spectra. 
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[Na/Fe] 

Fig. 12. [Ni/Fe] versus [Na/Fe] with the same symbols as used 
in Fig. |7| The line corresponds to the linear fit of [Ni/Fe] as a 
function of [Na/Fe] given in Eq. (2). 

stars with the lowest a/Fe ratios, which agrees with the trends 
seen for the solar twin stars. 

As discussed by Venn et al. (12004b . the Ni-Na correlation 
may be due to the fact that the yields of Na and the domi¬ 
nant isotope of Ni (Ni^^) produced in massive stars exploding 
as Type II SNe both depend on the neutron excess, which itself 
depends on metallicity and the a/Fe ratio. Neutrons are released 
by the ^^Ne(Qr,n)^^Mg reaction, where ^^Ne comes from double 
a-capture on ^"^N made from initial C and O in the CNO cycle. 
It is, however, difficult to understand why [Na/Fe] and [Ni/Fe] 
are not well correlated with stellar age like the other abundance 
ratios; it was argued above that a-capture elements produced in 
Type II SNe must be efficiently mixed in interstellar gas to ex¬ 
plain the small scatter of [Mg/Fe] at a given age, so why should 
there be variations of the neutron excess? 

5.2. [X/Fe]-7’c correlations 

As discussed in Sect. [1] sequestration of refractory elements in 
terrestrial planets or on interstellar dust particles may introduce 
star-to-star differences of abundance ratios, which are correlated 
with elemental condensation temperature Tq (Fodders 12003b . In 
order to see if this is the case, [X/Fe] has been plotted as a func¬ 
tion of Tc for all stars and linear-least squares fits with [X/Fe] 
weighted by the inverse square of its error have been obtained. 
The coefficients and the reduced of the fits are given in Table 
|7| Figure [13] shows some typical example^ 

As seen from Table]?] the reducedof the fits is larger than 
one for all stars indicating that there is never a perfect correla¬ 
tion between [X/Fe] and Tc. Nevertheless, for several stars the 
slope coefficient is more than three times larger than the error 
of the slope. The first star, HD 134664, in Fig. [13] with a slope 
A[X/Fe]/Arc = 6.7 ± 1.6 x lO^^dexK^' represents such a case. 
It has a ratio between refractory elements with Tq > 1200 K and 
the volatile elements C and O about 0.08 dex higher than in the 
Sun. Assuming that HD 134664 has a pristine composition, this 
means that the Sun has been depleted in refractory elements by 
as much as 20% as first advocated by Melendez et al. (12009b 
from abundance ratios in 11 solar twins. 

The second star in Fig. [13] HD 146233 (18 Sco) represents an 
intermediate case with A[X/Fe]/A7’c = 3.8+1.1 xlO^^dexK"^ 
This star has recently been studied in detail by Melendez et al. 

^ Yttrium is not included in the [X/Fel-Tc fits, because for many stars 
[Y/Fe] has a large deviation. 
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Table 7. Linear fits, [X/Fe] - a + b ■ Tq- 


Star 

a 

[dex] 

b 

[lO^MexK-*] 

2 

^re(d 

HD 2071 

-0.043 ±0.019 

3.1 ± 1.5 

4.9 

HD 8406 

-0.068 ±0.030 

4.4 ±2.3 

11.2 

HD 20782 

-0.060 ±0.034 

3.7 ±2.6 

16.3 

HD 27063 

-0.084 ±0.023 

5.5 ± 1.8 

6.8 

HD 28471 

-0.041 ±0.023 

4.3 ± 1.8 

6.8 

HD 38277 

-0.051 ±0.034 

2.9 ±2.6 

14.8 

HD 45184 

-0.069 ±0.013 

4.6 ± 1.0 

2.2 

HD 45289* 

+0.021 ±0.052 

1.9 ±4.0 

35.1 

HD 71334 

-0.059 ±0.026 

5.5 ±2.0 

8.3 

HD 78429 

-0.038 ±0.024 

2.8 ± 1.8 

8.4 

HD 88084 

-0.011 ±0.019 

1.9 ± 1.5 

5.0 

HD 92719 

-0.081 ±0.027 

5.4 ±2.1 

9.9 

HD 96116 

-0.113 ±0.042 

6.1 ±3.2 

21.7 

HD 96423 

-0.060 ±0.015 

5.5 ± 1.1 

3.1 

HD 134664 

-0.106 ±0.020 

6.7 ± 1.6 

5.6 

HD 146233 

-0.058 ±0.014 

3.8 ± 1.1 

2.7 

HD 183658 

+0.003 ±0.012 

0.4 ±0.9 

1.7 

HD 208704 

-0.041 ±0.029 

2.7 ±2.2 

11.5 

HD210918* 

+0.003 ±0.056 

1.4 ±4.3 

46.2 

HD 220507* 

+0.085 ±0.053 

-1.7 + 4.1 

43.9 

HD 222582 

-0.002 ±0.019 

1.7 ± 1.5 

5.4 


Notes. a-enhanced star. 


(I2014ab . who for 20 elements with Z < 30 find a Tc-slope of 
3.5 X lO^^dexK * in excellent agreement with the value found 
is this paperQ. Neutron capture elements were, however, found 
to have higher abundances than expected from the [X/FeJ-Tc re¬ 
lation. This is confirmed by deriving abundances of Y, Ba and 
Eu from the HARPS spectrum of HD 146233. 

The third star in Fig. [13] HD 183658, is the only one in the 
sample with a well-defined relation between [X/Fe] and Tq and a 
slope close to zero. Thus, it is the only star for which the [X/Fe] 
distribution matches that of the Sun, confirming that the Sun is 
unusual among solar twin stars. 

For the three cr-enhanced stars, there is no correlation be¬ 
tween [X/Fe] and Fc as shown in the case of HD 210918 in 
the last panel of Fig. [T3| In addition, there are hve stars with 
Yred >10 for which the correlation is also poor, mainly because 
[Na/Fe] and [Ni/Fe] deviate strongly from the mean trend. As 
shown in Fig.[T3]for the case of HD 20782 very different slopes 
are derived depending on whether one includes all elements or 
ht only elements with a condensation temperature Tq > 900 K as 
in Ramfrez et al. (120091) . Another solar twin (HIP 102152) with 
deviating [Na/Fe] and [Ni/Fe] values was analyzed by Monroe 
et al. (120131 ). who suggest that the low sodium (and nitrogen) 
abundance could be due to it’s high age (8 Gyr). However, stars 
with low Na and Ni abundances in this paper cover a broad 
range in age. Instead, it should be noted that HIP 102152 with 
[Na/Fe] = -0.044 and [Ni/Fe] = -0.026 according to Monroe 
et al. fits the Ni-Na relation shown in Fig. [T2| almost exactly. 

Recently, Adibekyan et al. (120141) have found evidence for 
a dependence of the [X/Fe]-Fc slope on stellar age using data 
based on HARPS spectra for 148 solar-like stars with 5600 < 
Feff < 6375 K, 4.10 < log g < 4.65 and -0.3 < [Fe/H] < h- 0.4. 
For this sample, the Fc-slope changes by about 6 x lO^^dex 


’ For 13 elements in common, the mean deviation and rms scatter 
for HD 146233 are A[X/H] (Melendez - this paper) = 0.003 ± 0.007 dex 
confirming the high precision estimated in the two papers. 






0 200 400 600 800 1000 1200 1400 1600 1800 

Tcond (K) 


Fig. 13. The relation between [X/Fe] and elemental condensa¬ 
tion temperature (Lodders 120031) for five representative stars. 
The lines show linear fits to the weighted values of [X/Fe]. In 
the case of HD 20782, the dashed line shows the fit for elements 
with Tc> 900 K. 


for an age change of 10 Gyr. A similar age dependence of the 
[X/Fe]-rc slope were found for a group of 59 metal-rich solar- 
analog stars by Ramfrez et al. (I2014al ). As seen from Fig. [T4| 
there is also evidence for an age dependence of the slopes de¬ 
rived in this paper. From a maximum likelihood fit to the data 
for the thin-disk stars I get 

b = 5.6 (+0.8) - 0.43 (±0.15) Age [Gyr], (4) 

where b is the [X/Fe] versus Tq slope in units of lO^^dexK''. 
The fit has a satisfactory= 1.5, but one star, HD 183658, and 
the Sun have a 3-sigma deviation. Within the errors, the change 
of the Fc-slope with age agrees with the result of Adibekyan 
et al.0. Furthermore, it is noted that the change in Fc-slope is 

* Note that Adibekyan et al. d2014> define the Fc-slope with the op¬ 
posite sign of that used in the present paper. 
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Fig. 14. The relation between the [X/FeJ-Tc slope and stellar age 
for thin-disk stars. The line shows the fit given in Eq. (3). 

mainly caused by the fact that the change of [C/Fe] and [O/Fe] 
with age is steeper than the corresponding change for refractory 
elements. 

6. Summary and conclusions 

In this paper, HARPS spectra with S /N > 600 have been used 
to derive very precise atmospheric parameters (cr(7’eff) ^ 6K; 
cr{logg) ^ 0.012) and elemental abundances for a sample of 
21 solar twin stars in the solar neighborhood. Small differen¬ 
tial non-FTE effects were taken into account when determining 
[Na/Fe] ratios and estimated to be negligible for the other abun¬ 
dance ratios derived. Differential 3D effects seem to be unim¬ 
portant, although the impact of possible variations in the mag¬ 
netic field strength deserves further studies. Altogether, the er¬ 
ror of the abundance ratios derived is estimated to be on the or¬ 
der of +0.01 dex. This high precision is supported by comparing 
abundances of HD 146233 (18 Sco) with those determined by 
Melendez et al. (I2014al l from UVES spectra. 

In addition, precise stellar ages are obtained by interpolat¬ 
ing between Yonsei-Yale isochrones in the logg-Tefr diagram. 
A comparison with ages derived by Ramirez et al. (I2014bl l in¬ 
dicates that the ages have internal errors ranging from 0.4 to 
0.8 Gyr. 

It is confirmed that the Sun has a lower ratio between re¬ 
fractory and volatile elements than the large majority of solar 
twin stars as found by Melendez et al. (I2009t . Only one star 
(HD 183658) has the same relation between [X/Fe] and elemen¬ 
tal condensation temperature Tc as the Sun. For several of the 
other stars, there is a positive correlation between [X/Fe] and Tq 
with slope coefficients ranging by more than a factor of two. This 
may be explained by various degrees of sequestration of refrac¬ 
tory elements on terrestrial planets (Melendez et al. 120091 1 or by 
dust-gas segregation in circumstellar disks (Gaidos l2015ll . 

For most of the solar twins, the correlation between [X/Fe] 
and Tq is rather poor, i.e., > 1 for the linear fits. For many 

of the elements there is, on the other hand, an astonishingly tight 
correlation between [X/Fe] and stellar age with amplitudes up 
to ~ 0.15 dex over an age interval of 8 Gyr. These age correla¬ 
tions may be related to an increasing number of Type la SNe 
relative to the number of Type II SNe as time goes on and/or 
an evolving initial mass function in the Galactic disk. If so, the 
slopes of of the age relations for the various elements provide 
new constraints on supernovae yields and Galactic chemical evo¬ 
lution. Furthermore, the small scatter of [X/Fe] in the age re¬ 
lations suggests that nucleosynthesis products of SNe are well 


mixed locally before new generations of stars are formed. In 
contrast, there is a large scatter of [Fe/H] at a given age, which 
may be explained by infall of metal-poor gas clouds followed 
by star formation before mixing evens out the inhomogeneities 
(Edvardsson et al. ll993l l. 

[Ca/Fe] is nearly constant as a function of stellar age in con¬ 
trast to an increasing [X/Fe] for the other a-capture elements. 
This cannot be explained by yields calculated for supernovae 
of Types II and la, and perhaps it indicates that the new class 
of low-luminosity Ca-rich supernovae discussed by Perets et al. 
(120101) play an important role for the chemical evolution of the 
Galactic disk. 

[Y/Mg] appears to be a sensitive chronometer for Galactic 
evolution as seen from Fig. [TOl If [Y/Mg] can be measured with 
a precision of 0.04 dex, the age of solar metallicity stars can be 
estimated from Eq. (2) with a precision of 1 Gyr, but further stud¬ 
ies are needed to see how this relation depends on [Fe/H]. 

The C/O ratio is found to evolve very little over the life 
time of the Galactic disk, i.e., to lie within +0.04 from the so¬ 
lar value of C/O ^ 0.55; in contrast [C/Fe] and [O/Fe] change by 
~ 0.15 dex. This is of particular interest for discussions of com¬ 
positions of exoplanets, because C/O is a key determinant for 
the composition of solids that condense from gas in circumstel¬ 
lar disks. 

Na and Ni are exceptions to the tight correlation between 
[X/Fe] and stellar age, but interestingly there is a strong correla¬ 
tion between [Ni/Fe] and [Na/Fe] with a slope of ^ 0.6 and an 
amplitude of ~ 0.07 dex in [Ni/Fe]. A similar Ni-Na relation has 
previously been found for more metal-poor halo stars (Nissen & 
Schuster 1 1 9971 l20T0l l and has been ascribed to the fact that the 
yields of Na and Ni produced in massive stars and expelled by 
Type II SNe both depend on the neutron excess. It is, however, 
difficult to understand why [Na/Fe] and [Ni/Fe] are not well cor¬ 
related with stellar age like the other abundance ratios. 

The changes of [X/Fe] with stellar age and the correlated 
variations of [Ni/Fe] and [Na/Fe] among solar twin stars com¬ 
plicate the use of the [X/Fe]-7’c slope as a possible signature of 
terrestrial planets around stars. Altogether, the interpretation of 
abundance ratios in solar twin stars stars, when determined with 
a precision of ~ 0.01 dex, is difficult, because [X/Fe] seems to 
be affected both by planet formation or dust-gas segregation and 
chemical evolution in the Galactic disk. Further studies of high- 
precision differential abundances for solar twins and stars with 
other metallicities including wide binaries are needed. 
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Table 2. List of spectral lines. 


Element 

Wavelength 

[A] 

Exc. pot. 
[eV] 

EWo 

[mA] 

Ci 

5052.15 

7.685 

36.0 

Ci 

5380.32 

7.685 

21.4 

[Oi] 

6300.31 

0.000 

3.6 

Nai 

6154.23 

2.102 

38.2 

Nai 

6160.75 

2.104 

58.8 

Mg I 

4730.04 

4.340 

71.5 

Mg I 

5711.10 

4.345 

105.6 

Ah 

6696.03 

3.143 

37.7 

Ah 

6698.67 

3.143 

21.0 

Si I 

5517.55 

5.080 

13.8 

Si I 

5645.62 

4.929 

36.6 

Si I 

5665.56 

4.920 

41.2 

Si I 

5793.08 

4.929 

43.9 

Si I 

6125.03 

5.614 

32.5 

Si I 

6145.02 

5.616 

38.8 

Si I 

6243.82 

5.616 

48.1 

Si I 

6244.48 

5.616 

46.3 

Si I 

6721.84 

5.862 

44.7 

Si I 

6741.63 

5.984 

15.9 

Si 

6045.98 

7.868 

18.3 

Si 

6052.67 

7.870 

11.8 

Si 

6743.58 

7.866 

8.6 

Si 

6757.14 

7.870 

18.7 

Cai 

5260.39 

2.521 

32.8 

Cai 

5512.99 

2.933 

88.3 

Cai 

5581.98 

2.523 

96.5 

Cai 

5590.13 

2.521 

92.6 

Cai 

5867.57 

2.933 

25.0 

Cai 

6166.44 

2.521 

71.1 

Cai 

6455.60 

2.523 

57.4 

Til 

4913.62 

1.873 

51.5 

Til 

5113.45 

1.443 

27.9 

Til 

5219.71 

0.021 

28.1 

Til 

5295.78 

1.067 

12.7 

Til 

5490.16 

1.460 

22.2 

Til 

5739.48 

2.249 

8.0 

Til 

5866.46 

1.066 

48.6 

Til 

6091.18 

2.267 

15.2 

Til 

6126.22 

1.066 

22.6 

Til 

6258.11 

1.443 

52.2 

Til 

6261.11 

1.429 

48.5 

Till 

5211.54 

2.590 

33.7 

Till 

5381.03 

1.565 

60.9 

Till 

5418.77 

1.582 

49.0 

Cri 

5214.13 

3.369 

17.0 

Cri 

5238.97 

2.709 

16.9 

Cri 

5247.57 

0.960 

83.4 

Cri 

5272.00 

3.449 

23.8 

Cri 

5287.18 

3.438 

11.3 

Cri 

5296.70 

0.983 

93.9 

Cri 

5348.33 

1.004 

100.2 

Cri 

5783.07 

3.323 

31.4 

Cri 

5783.87 

3.322 

44.6 

Cri 

6661.08 

4.193 

13.0 

Cm 

5237.33 

4.073 

53.5 

Cm 

5246.78 

3.714 

16.2 

Fei 

5295.32 

4.415 

29.9 

Fei 

5373.71 

4.473 

63.9 

Fei 

5379.58 

3.694 

61.9 
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Table 2. continued 


Element 

Wavelength 

(A) 

Exc. pot. 
(eV) 

FWo 

(mA) 

Fei 

5386.34 

4.154 

33.1 

Fei 

5466.99 

3.573 

35.9 

Fei 

5522.45 

4.209 

44.2 

Fei 

5546.51 

4.371 

52.9 

Fei 

5560.22 

4.434 

52.6 

Fei 

5577.03 

5.033 

11.6 

Fei 

5618.64 

4.209 

50.4 

Fei 

5636.71 

3.640 

19.6 

Fei 

5650.00 

5.100 

36.5 

Fei 

5651.48 

4.473 

19.0 

Fei 

5661.35 

4.284 

23.0 

Fei 

5679.03 

4.652 

60.4 

Fei 

5705.47 

4.301 

38.4 

Fei 

5855.09 

4.608 

22.7 

Fei 

6079.02 

4.652 

46.5 

Fei 

6082.72 

2.223 

35.4 

Fei 

6093.65 

4.607 

31.1 

Fei 

6096.67 

3.984 

37.8 

Fei 

6151.62 

2.176 

50.1 

Fei 

6157.73 

4.076 

62.4 

Fei 

6165.36 

4.143 

45.0 

Fei 

6173.34 

2.223 

68.7 

Fei 

6188.00 

3.943 

48.4 

Fei 

6226.74 

3.883 

29.6 

Fei 

6240.65 

2.223 

49.0 

Fei 

6270.23 

2.858 

52.2 

Fei 

6271.28 

3.332 

24.4 

Fei 

6380.75 

4.186 

52.7 

Fei 

6392.54 

2.279 

17.6 

Fei 

6498.94 

0.958 

47.1 

Fei 

6597.57 

4.795 

44.0 

Fei 

6625.04 

1.011 

15.6 

Fei 

6703.58 

2.759 

37.3 

Fei 

6705.10 

4.607 

47.6 

Fei 

6710.32 

1.485 

16.0 

Fei 

6713.75 

4.795 

21.3 

Fei 

6725.36 

4.103 

17.8 

Fei 

6726.67 

4.607 

47.3 

Fei 

6733.15 

4.638 

26.9 

Fei 

6739.52 

1.557 

12.0 

Fei 

6793.27 

4.076 

13.0 

Fei 

6806.86 

2.727 

34.9 

Fei 

6810.27 

4.607 

50.1 

Fei 

6843.65 

4.548 

61.5 

Fell 

5414.08 

3.222 

27.3 

Fell 

5425.26 

3.200 

41.5 

Fell 

6084.10 

3.200 

21.2 

Fell 

6149.25 

3.889 

36.7 

Fell 

6247.56 

3.892 

53.1 

Fell 

6369.46 

2.891 

18.8 

Fell 

6416.93 

3.892 

40.1 

Fell 

6432.68 

2.892 

41.6 

Fell 

6456.39 

3.904 

63.4 

Nil 

4953.21 

3.740 

56.2 

Nil 

5010.94 

3.635 

49.9 

Nil 

5643.09 

4.164 

15.9 

Nil 

5805.23 

4.167 

41.5 

Nil 

6086.29 

4.266 

44.2 

Nil 

6108.12 

1.676 

65.6 
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Table 2. continued 


Element 

Wavelength 

(A) 

Exc. pot. 
(eV) 

EWo 

(mA) 

Nil 

6130.14 

4.266 

22.3 

Nil 

6176.82 

4.088 

64.1 

Nil 

6177.25 

1.826 

14.8 

Nil 

6204.61 

4.088 

22.1 

Nil 

6378.26 

4.154 

32.8 

Nil 

6643.64 

1.676 

94.7 

Nil 

6767.78 

1.826 

79.7 

Nil 

6772.32 

3.657 

50.1 

Zni 

4722.16 

4.030 

71.6 

Zni 

4810.54 

4.080 

75.0 

Zni 

6362.35 

5.790 

20.4 

Yii 

4883.69 

1.084 

59.5 

Yii 

5087.43 

1.084 

48.7 

Yii 

5200.42 

0.992 

38.1 
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Table 3. Atmospheric parameters and abundance ratios. 


Star Teff log 9 [Fe/H] ^turb [C/Fe] [0/Fe] [Na/Fe] [Mg/Fe] [Al/Fe] [Si/Fe] [S/Fe] [Ca/Fe] [Ti/Fe] [Cr/Fe] [Ni/Fe] [Zn/Fe] [Y/Fe] 

K km s“' 


HD2071 5724 4.490 -0.084 0.96 -0.022 -0.022 -0.032 0.001 

HD8406 5730 4.479 -0.105 0.95 -0.036 -0.028 -0.059 0.000 

HD20782 5776 4.345 -0.058 1.04 -0.019 0.024 -0.086 0.011 

HD27063 5779 4.469 0.064 0.99 -0.082 -0.078 -0.041 -0.018 

HD28471 5754 4.380 -0.054 1.02 -0.023 0.006 -0.013 0.037 

HD38277 5860 4.270 -0.070 1.17 -0.012 0.049 -0.078 0.012 

HD45184 5871 4.445 0.047 1.06 -0.069 -0.055 -0.031 -0.009 

HD45289* 5718 4.284 -0.020 1.06 0.058 0.146 -0.009 0.100 

HD71334 5701 4.374 -0.075 0.98 -0.024 0.031 -0.034 0.039 

HD78429 5756 4.272 0.078 1.05 -0.025 0.004 -0.055 0.022 

HD88084 5768 4.424 -0.091 1.02 0.017 -0.010 0.003 0.028 

HD92719 5813 4.488 -0.112 1.00 -0.060 -0.009 -0.050 0.002 

HD96116 5846 4.503 0.006 1.02 -0.088 -0.046 -0.100 -0.041 

HD96423 5714 4.359 0.113 0.99 -0.056 -0.076 0.008 0.021 

HD134664 5853 4.452 0.093 1.01 -0.116 -0.087 -0.058 -0.022 

HD146233 5809 4.434 0.046 1.02 -0.066 -0.028 -0.028 -0.010 

HD183658 5809 4.402 0.035 1.04 -0.002 -0.040 0.024 0.009 

HD208704 5828 4.346 -0.091 1.08 -0.009 0.033 -0.062 0.017 

HD210918* 5748 4.319 -0.095 1.07 0.054 0.146 -0.057 0.072 

HD220507* 5690 4.247 0.013 1.07 0.102 0.159 0.036 0.126 

HD222582 5784 4.361 -0.004 1.07 -0.001 0.043 0.014 0.037 


-0.002 0.001 -0.007 0.024 0.019 0.005 -0.029 -0.026 0.040 

-0.020 -0.010 -0.021 0.031 0.009 0.010 -0.041 -0.044 0.050 

0.015 -0.006 -0.013 0.015 0.028 -0.004 -0.042 -0.019 -0.110 

-0.029 -0.019 -0.039 0.023 -0.005 0.016 -0.022 -0.055 0.067 

0.065 0.022 -0.010 0.026 0.046 -0.007 -0.005 0.018 -0.067 

0.012 0.006 -0.017 0.007 0.023 -0.018 -0.042 -0.016 -0.129 

-0.015 -0.005 -0.045 0.015 0.011 -0.002 -0.020 -0.030 0.040 

0.128 0.055 0.056 0.049 0.100 0.002 0.006 0.096 -0.094 

0.063 0.019 -0.003 0.034 0.041 0.000 -0.014 -0.007 -0.076 

0.034 -0.002 0.001 0.006 0.019 0.005 -0.020 0.014 -0.101 

0.053 0.021 0.008 0.019 0.035 -0.009 -0.003 0.009 -0.044 

-0.014 0.003 -0.062 0.031 0.021 -0.001 -0.044 -0.046 0.042 

-0.063 -0.036 -0.052 0.019 -0.019 0.013 -0.056 -0.085 0.095 

0.056 0.005 -0.035 0.005 0.028 0.005 0.021 -0.014 -0.033 

-0.035 -0.019 -0.048 0.009 0.004 0.006 -0.022 -0.055 0.050 

-0.020 -0.007 -0.029 0.010 0.009 0.010 -0.013 -0.036 0.041 

0.018 0.009 0.006 0.002 0.008 -0.005 0.013 0.008 -0.025 

0.023 0.005 -0.020 0.015 0.024 -0.012 -0.033 0.000 -0.111 

0.102 0.029 0.032 0.043 0.080 -0.012 -0.029 0.034 -0.077 

0.153 0.071 0.084 0.051 0.114 0.004 0.017 0.121 -0.085 

0.058 0.028 -0.021 0.015 0.041 0.000 0.005 0.025 -0.065 


Notes. *** a-enhanced star. 







